A view of the pan-genome of domesticated Cowpea (Vigna unguiculata [L.] Walp.).

Cowpea, Vigna unguiculata L. Walp., is a diploid warm-season legume of critical importance as both food and fodder in sub-Saharan Africa. This species is also grown in Northern Africa, Europe, Latin America, North America, and East to Southeast Asia. To capture the genomic diversity of domesticates of this important legume, de novo genome assemblies were produced for representatives of six subpopulations of cultivated cowpea identified previously from genotyping of several hundred diverse accessions. In the most complete assembly (IT97K-499-35), 26,026 core and 4963 noncore genes were identified, with 35,436 pan genes when considering all seven accessions. GO terms associated with response to stress and defense response were highly enriched among the noncore genes, while core genes were enriched in terms related to transcription factor activity, and transport and metabolic processes. Over 5 million single nucleotide polymorphisms (SNPs) relative to each assembly and over 40 structural variants >1 Mb in size were identified by comparing genomes. Vu10 was the chromosome with the highest frequency of SNPs, and Vu04 had the most structural variants. Noncore genes harbor a larger proportion of potentially disruptive variants than core genes, including missense, stop gain, and frameshift mutations; this suggests that noncore genes substantially contribute to diversity within domesticated cowpea.

UDP-glucosyltransferase HvUGT13248 confers type II resistance to Fusarium graminearum in barley.

Fusarium head blight (FHB) of barley (Hordeum vulgare) causes yield losses and accumulation of trichothecene mycotoxins (e.g. deoxynivalenol [DON]) in grains. Glucosylation of DON to the nontoxic DON-3-O-glucoside (D3G) is catalyzed by UDP-glucosyltransferases (UGTs), such as barley UGT13248. We explored the natural diversity of UGT13248 in 496 barley accessions and showed that all carried potential functional alleles of UGT13248, as no genotypes showed strongly increased seedling sensitivity to DON. From a TILLING population, we identified 2 mutant alleles (T368I and H369Y) that, based on protein modeling, likely affect the UDP-glucose binding of UGT13248. In DON feeding experiments, DON-to-D3G conversion was strongly reduced in spikes of these mutants compared to controls, and plants overexpressing UGT13248 showed increased resistance to DON and increased DON-to-D3G conversion. Moreover, field-grown plants carrying the T368I or H369Y mutations inoculated with Fusarium graminearum showed increased FHB disease severity and reduced D3G production. Barley is generally considered to have type II resistance that limits the spread of F. graminearum from the infected spikelet to adjacent spikelets. Point inoculation experiments with F. graminearum showed increased infection spread in T368I and H369Y across the spike compared to wild type, while overexpression plants showed decreased spread of FHB symptoms. Confocal microscopy revealed that F. graminearum spread to distant rachis nodes in T368I and H369Y mutants but was arrested at the rachis node of the inoculated spikelet in wild-type plants. Taken together, our data reveal that UGT13248 confers type II resistance to FHB in barley via conjugation of DON to D3G.

Diamonds in the not-so-rough: Wild relative diversity hidden in crop genomes.

Crop production is becoming an increasing challenge as the global population grows and the climate changes. Modern cultivated crop species are selected for productivity under optimal growth environments and have often lost genetic variants that could allow them to adapt to diverse, and now rapidly changing, environments. These genetic variants are often present in their closest wild relatives, but so are less desirable traits. How to preserve and effectively utilize the rich genetic resources that crop wild relatives offer while avoiding detrimental variants and maladaptive genetic contributions is a central challenge for ongoing crop improvement. This Essay explores this challenge and potential paths that could lead to a solution.

Unloading potatoes: Potato breeding moves forward with only half the genome.

Potato is one of the most important food crops in the world and, in contrast to other staples, has not seen large improvements in yield. Agha, Shannon, and Morrell preview an article recently published in Cell, "Phylogenomic discovery of deleterious mutations facilitates hybrid potato breeding," which advances potato breeding strategies via a genetic approach.

Genomic conservation of crop wild relatives: A case study of citrus.

Conservation of crop wild relatives is critical for plant breeding and food security. The lack of clarity on the genetic factors that lead to endangered status or extinction create difficulties when attempting to develop concrete recommendations for conserving a citrus wild relative: the wild relatives of crops. Here, we evaluate the conservation of wild kumquat (Fortunella hindsii) using genomic, geographical, environmental, and phenotypic data, and forward simulations. Genome resequencing data from 73 accessions from the Fortunella genus were combined to investigate population structure, demography, inbreeding, introgression, and genetic load. Population structure was correlated with reproductive type (i.e., sexual and apomictic) and with a significant differentiation within the sexually reproducing population. The effective population size for one of the sexually reproducing subpopulations has recently declined to ~1,000, resulting in high levels of inbreeding. In particular, we found that 58% of the ecological niche overlapped between wild and cultivated populations and that there was extensive introgression into wild samples from cultivated populations. Interestingly, the introgression pattern and accumulation of genetic load may be influenced by the type of reproduction. In wild apomictic samples, the introgressed regions were primarily heterozygous, and genome-wide deleterious variants were hidden in the heterozygous state. In contrast, wild sexually reproducing samples carried a higher recessive deleterious burden. Furthermore, we also found that sexually reproducing samples were self-incompatible, which prevented the reduction of genetic diversity by selfing. Our population genomic analyses provide specific recommendations for distinct reproductive types and monitoring during conservation. This study highlights the genomic landscape of a wild relative of citrus and provides recommendations for the conservation of crop wild relatives.

Selection and adaptive introgression guided the complex evolutionary history of the European common bean.

Domesticated crops have been disseminated by humans over vast geographic areas. Common bean (Phaseolus vulgaris L.) was introduced in Europe after 1492. Here, by combining whole-genome profiling, metabolic fingerprinting and phenotypic characterisation, we show that the first common bean cultigens successfully introduced into Europe were of Andean origin, after Francisco Pizarro's expedition to northern Peru in 1529. We reveal that hybridisation, selection and recombination have shaped the genomic diversity of the European common bean in parallel with political constraints. There is clear evidence of adaptive introgression into the Mesoamerican-derived European genotypes, with 44 Andean introgressed genomic segments shared by more than 90% of European accessions and distributed across all chromosomes except PvChr11. Genomic scans for signatures of selection highlight the role of genes relevant to flowering and environmental adaptation, suggesting that introgression has been crucial for the dissemination of this tropical crop to the temperate regions of Europe.

The evolutionary patterns of barley pericentromeric chromosome regions, as shaped by linkage disequilibrium and domestication.

The distribution of recombination events along large cereal chromosomes is uneven and is generally restricted to gene-rich telomeric ends. To understand how the lack of recombination affects diversity in the large pericentromeric regions, we analysed deep exome capture data from a final panel of 815 Hordeum vulgare (barley) cultivars, landraces and wild barleys, sampled from across their eco-geographical ranges. We defined and compared variant data across the pericentromeric and non-pericentromeric regions, observing a clear partitioning of diversity both within and between chromosomes and germplasm groups. Dramatically reduced diversity was found in the pericentromeres of both cultivars and landraces when compared with wild barley. We observed a mixture of completely and partially differentiated single-nucleotide polymorphisms (SNPs) between domesticated and wild gene pools, suggesting that domesticated gene pools were derived from multiple wild ancestors. Patterns of genome-wide linkage disequilibrium, haplotype block size and number, and variant frequency within blocks showed clear contrasts among individual chromosomes and between cultivars and wild barleys. Although most cultivar chromosomes shared a single major pericentromeric haplotype, chromosome 7H clearly differentiated the two-row and six-row types associated with different geographical origins. Within the pericentromeric regions we identified 22 387 non-synonymous SNPs, 92 of which were fixed for alternative alleles in cultivar versus wild accessions. Surprisingly, only 29 SNPs found exclusively in the cultivars were predicted to be 'highly deleterious'. Overall, our data reveal an unconventional pericentromeric genetic landscape among distinct barley gene pools, with different evolutionary processes driving domestication and diversification.

Genomic variants affecting homoeologous gene expression dosage contribute to agronomic trait variation in allopolyploid wheat.

Allopolyploidy greatly expands the range of possible regulatory interactions among functionally redundant homoeologous genes. However, connection between the emerging regulatory complexity and expression and phenotypic diversity in polyploid crops remains elusive. Here, we use diverse wheat accessions to map expression quantitative trait loci (eQTL) and evaluate their effects on the population-scale variation in homoeolog expression dosage. The relative contribution of cis- and trans-eQTL to homoeolog expression variation is strongly affected by both selection and demographic events. Though trans-acting effects play major role in expression regulation, the expression dosage of homoeologs is largely influenced by cis-acting variants, which appear to be subjected to selection. The frequency and expression of homoeologous gene alleles showing strong expression dosage bias are predictive of variation in yield-related traits, and have likely been impacted by breeding for increased productivity. Our study highlights the importance of genomic variants affecting homoeolog expression dosage in shaping agronomic phenotypes and points at their potential utility for improving yield in polyploid crops.

Fast neutron mutagenesis in soybean enriches for small indels and creates frameshift mutations.

The mutagenic effects of ionizing radiation have been used for decades to create novel variants in experimental populations. Fast neutron (FN) bombardment as a mutagen has been especially widespread in plants, with extensive reports describing the induction of large structural variants, i.e., deletions, insertions, inversions, and translocations. However, the full spectrum of FN-induced mutations is poorly understood. We contrast small insertions and deletions (indels) observed in 27 soybean lines subject to FN irradiation with the standing indels identified in 107 diverse soybean lines. We use the same populations to contrast the nature and context (bases flanking a nucleotide change) of single-nucleotide variants. The accumulation of new single-nucleotide changes in FN lines is marginally higher than expected based on spontaneous mutation. In FN-treated lines and in standing variation, C→T transitions and the corresponding reverse complement G→A transitions are the most abundant and occur most frequently in a CpG local context. These data indicate that most SNPs identified in FN lines are likely derived from spontaneous de novo processes in generations following mutagenesis rather than from the FN irradiation mutagen. However, small indels in FN lines differ from standing variants. Short insertions, from 1 to 6 bp, are less abundant than in standing variation. Short deletions are more abundant and prone to induce frameshift mutations that should disrupt the structure and function of encoded proteins. These findings indicate that FN irradiation generates numerous small indels, increasing the abundance of loss-of-function mutations that impact single genes.

Plant Genome Editing and the Relevance of Off-Target Changes.

Site-directed nucleases (SDNs) used for targeted genome editing are powerful new tools to introduce precise genetic changes into plants. Like traditional approaches, such as conventional crossing and induced mutagenesis, genome editing aims to improve crop yield and nutrition. Next-generation sequencing studies demonstrate that across their genomes, populations of crop species typically carry millions of single nucleotide polymorphisms and many copy number and structural variants. Spontaneous mutations occur at rates of ∼10-8 to 10-9 per site per generation, while variation induced by chemical treatment or ionizing radiation results in higher mutation rates. In the context of SDNs, an off-target change or edit is an unintended, nonspecific mutation occurring at a site with sequence similarity to the targeted edit region. SDN-mediated off-target changes can contribute to a small number of additional genetic variants compared to those that occur naturally in breeding populations or are introduced by induced-mutagenesis methods. Recent studies show that using computational algorithms to design genome editing reagents can mitigate off-target edits in plants. Finally, crops are subject to strong selection to eliminate off-type plants through well-established multigenerational breeding, selection, and commercial variety development practices. Within this context, off-target edits in crops present no new safety concerns compared to other breeding practices. The current generation of genome editing technologies is already proving useful to develop new plant varieties with consumer and farmer benefits. Genome editing will likely undergo improved editing specificity along with new developments in SDN delivery and increasing genomic characterization, further improving reagent design and application.

The Fate of Deleterious Variants in a Barley Genomic Prediction Population.

Targeted identification and purging of deleterious genetic variants has been proposed as a novel approach to animal and plant breeding. This strategy is motivated, in part, by the observation that demographic events and strong selection associated with cultivated species pose a "cost of domestication." This includes an increase in the proportion of genetic variants that are likely to reduce fitness. Recent advances in DNA resequencing and sequence constraint-based approaches to predict the functional impact of a mutation permit the identification of putatively deleterious SNPs (dSNPs) on a genome-wide scale. Using exome capture resequencing of 21 barley lines, we identified 3855 dSNPs among 497,754 total SNPs. We generated whole-genome resequencing data of Hordeum murinum ssp. glaucum as a phylogenetic outgroup to polarize SNPs as ancestral vs. derived. We also observed a higher proportion of dSNPs per synonymous SNPs (sSNPs) in low-recombination regions of the genome. Using 5215 progeny from a genomic prediction experiment, we examined the fate of dSNPs over three breeding cycles. Adjusting for initial frequency, derived alleles at dSNPs reduced in frequency or were lost more often than other classes of SNPs. The highest-yielding lines in the experiment, as chosen by standard genomic prediction approaches, carried fewer homozygous dSNPs than randomly sampled lines from the same progeny cycle. In the final cycle of the experiment, progeny selected by genomic prediction had a mean of 5.6% fewer homozygous dSNPs relative to randomly chosen progeny from the same cycle.

Environmental Association Identifies Candidates for Tolerance to Low Temperature and Drought.

Barley (Hordeum vulgare ssp. vulgare) is cultivated from the equator to the Arctic Circle. The wild progenitor species, Hordeum vulgare ssp. spontaneum, occupies a relatively narrow latitudinal range (∼30 - 40° N) primarily at low elevation (< 1,500 m). Adaptation to the range of cultivation has occurred over ∼8,000 years. The genetic basis of adaptation is amenable to study through environmental association. An advantage of environmental association in a well-characterized crop is that many loci that contribute to climatic adaptation and abiotic stress tolerance have already been identified. This provides the opportunity to determine if environmental association approaches effectively identify these loci of large effect. Using published genotyping from 7,864 SNPs in 803 barley landraces, we examined allele frequency differentiation across multiple partitions of the data and mixed model associations relative to bioclimatic variables. Using newly generated resequencing data from a subset of these landraces, we tested for linkage disequilibrium (LD) between SNPs queried in genotyping and SNPs in neighboring loci. Six loci previously reported to contribute to adaptive differences in flowering time and abiotic stress in barley and six loci previously identified in other plant species were identified in our analyses. In many cases, patterns of LD are consistent with the causative variant occurring in the immediate vicinity of the queried SNP. The identification of barley orthologs to well-characterized genes may provide a new understanding of the nature of adaptive variation and could permit a more targeted use of potentially adaptive variants in barley breeding and germplasm improvement.

Development of a Multiparent Population for Genetic Mapping and Allele Discovery in Six-Row Barley.

Germplasm collections hold valuable allelic diversity for crop improvement and genetic mapping of complex traits. To gain access to the genetic diversity within the USDA National Small Grain Collection (NSGC), we developed the Barley Recombinant Inbred Diverse Germplasm Population (BRIDG6), a six-row spring barley multiparent population (MPP) with 88 cultivated accessions crossed to a common parent (Rasmusson). The parents were randomly selected from a core subset of the NSGC that represents the genetic diversity of landrace and breeding accessions. In total, we generated 6160 F5 recombinant inbred lines (RILs), with an average of 69 and a range of 37-168 RILs per family, that were genotyped with 7773 SNPs, with an average of 3889 SNPs segregating per family. We detected 23 quantitative trait loci (QTL) associated with flowering time with five QTL found coincident with previously described flowering time genes. A major QTL was detected near the flowering time gene, HvPpd-H1 which affects photoperiod. Haplotype-based analysis of HvPpd-H1 identified private alleles to families of Asian origin conferring both positive and negative effects, providing the first observation of flowering time-related alleles private to Asian accessions. We evaluated several subsampling strategies to determine the effect of sample size on the power of QTL detection, and found that, for flowering time in barley, a sample size >50 families or 3000 individuals results in the highest power for QTL detection. This MPP will be useful for uncovering large and small effect QTL for traits of interest, and identifying and utilizing valuable alleles from the NSGC for barley improvement.

Publisher Correction: Exome sequencing highlights the role of wild-relative introgression in shaping the adaptive landscape of the wheat genome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Exome sequencing highlights the role of wild-relative introgression in shaping the adaptive landscape of the wheat genome.

Introgression is a potential source of beneficial genetic diversity. The contribution of introgression to adaptive evolution and improvement of wheat as it was disseminated worldwide remains unknown. We used targeted re-sequencing of 890 diverse accessions of hexaploid and tetraploid wheat to identify wild-relative introgression. Introgression, and selection for improvement and environmental adaptation, each reduced deleterious allele burden. Introgression increased diversity genome wide and in regions harboring major agronomic genes, and contributed alleles explaining a substantial proportion of phenotypic variation. These results suggest that historic gene flow from wild relatives made a substantial contribution to the adaptive diversity of modern bread wheat.

Comparative Genomics Approaches Accurately Predict Deleterious Variants in Plants.

Recent advances in genome resequencing have led to increased interest in prediction of the functional consequences of genetic variants. Variants at phylogenetically conserved sites are of particular interest, because they are more likely than variants at phylogenetically variable sites to have deleterious effects on fitness and contribute to phenotypic variation. Numerous comparative genomic approaches have been developed to predict deleterious variants, but the approaches are nearly always assessed based on their ability to identify known disease-causing mutations in humans. Determining the accuracy of deleterious variant predictions in nonhuman species is important to understanding evolution, domestication, and potentially to improving crop quality and yield. To examine our ability to predict deleterious variants in plants we generated a curated database of 2,910 Arabidopsis thaliana mutants with known phenotypes. We evaluated seven approaches and found that while all performed well, their relative ranking differed from prior benchmarks in humans. We conclude that deleterious mutations can be reliably predicted in A. thaliana and likely other plant species, but that the relative performance of various approaches does not necessarily translate from one species to another.

Genetic Costs of Domestication and Improvement.

The "cost of domestication" hypothesis posits that the process of domesticating wild species can result in an increase in the number, frequency, and/or proportion of deleterious genetic variants that are fixed or segregating in the genomes of domesticated species. This cost may limit the efficacy of selection and thus reduce genetic gains in breeding programs for these species. Understanding when and how deleterious mutations accumulate can also provide insight into fundamental questions about the interplay of demography and selection. Here we describe the evolutionary processes that may contribute to deleterious variation accrued during domestication and improvement, and review the available evidence for "the cost of domestication" in animal and plant genomes. We identify gaps and explore opportunities in this emerging field, and finally offer suggestions for researchers and breeders interested in understanding or avoiding the consequences of an increased number or frequency of deleterious variants in domesticated species.

Deleterious Variants in Asian Rice and the Potential Cost of Domestication.

Many SNPs are predicted to encode deleterious amino acid variants. These slightly deleterious mutations can provide unique insights into population history, the dynamics of selection, and the genetic bases of phenotypes. This is especially true for domesticated species, where a history of bottlenecks and selection may affect the frequency of deleterious variants and signal a "cost of domestication". Here, we investigated the numbers and frequencies of deleterious variants in Asian rice (Oryza sativa), focusing on two varieties (japonica and indica) and their wild relative (O. rufipogon). We investigated three signals of a potential cost of domestication in Asian rice relative to O. rufipogon: an increase in the frequency of deleterious SNPs (dSNPs), an enrichment of dSNPs compared with synonymous SNPs (sSNPs), and an increased number of deleterious variants. We found evidence for all three signals, and domesticated individuals contained ∼3-4% more deleterious alleles than wild individuals. Deleterious variants were enriched within low recombination regions of the genome and experienced frequency increases similar to sSNPs within regions of putative selective sweeps. A characteristic feature of rice domestication was a shift in mating system from outcrossing to predominantly selfing. Forward simulations suggest that this shift in mating system may have been the dominant factor in shaping both deleterious and neutral diversity in rice.

angsd-wrapper: utilities for analysing next-generation sequencing data.

High-throughput sequencing has changed many aspects of population genetics, molecular ecology and related fields, affecting both experimental design and data analysis. The software package angsd allows users to perform a number of population genetic analyses on high-throughput sequencing data. angsd uses probabilistic approaches which can directly make use of genotype likelihoods; thus, SNP calling is not required for comparative analyses. This takes advantage of all the sequencing data and produces more accurate results for samples with low sequencing depth. Here, we present angsd-wrapper, a set of wrapper scripts that provides a user-friendly interface for running angsd and visualizing results. angsd-wrapper supports multiple types of analyses including estimates of nucleotide sequence diversity neutrality tests, principal component analysis, estimation of admixture proportions for individual samples and calculation of statistics that quantify recent introgression. angsd-wrapper also provides interactive graphing of angsd results to enhance data exploration. We demonstrate the usefulness of angsd-wrapper by analysing resequencing data from populations of wild and domesticated Zea. angsd-wrapper is freely available from https://github.com/mojaveazure/angsd-wrapper.

The Role of Deleterious Substitutions in Crop Genomes.

Populations continually incur new mutations with fitness effects ranging from lethal to adaptive. While the distribution of fitness effects of new mutations is not directly observable, many mutations likely either have no effect on organismal fitness or are deleterious. Historically, it has been hypothesized that a population may carry many mildly deleterious variants as segregating variation, which reduces the mean absolute fitness of the population. Recent advances in sequencing technology and sequence conservation-based metrics for inferring the functional effect of a variant permit examination of the persistence of deleterious variants in populations. The issue of segregating deleterious variation is particularly important for crop improvement, because the demographic history of domestication and breeding allows deleterious variants to persist and reach moderate frequency, potentially reducing crop productivity. In this study, we use exome resequencing of 15 barley accessions and genome resequencing of 8 soybean accessions to investigate the prevalence of deleterious single nucleotide polymorphisms (SNPs) in the protein-coding regions of the genomes of two crops. We conclude that individual cultivars carry hundreds of deleterious SNPs on average, and that nonsense variants make up a minority of deleterious SNPs. Our approach identifies known phenotype-altering variants as deleterious more frequently than the genome-wide average, suggesting that putatively deleterious variants are likely to affect phenotypic variation. We also report the implementation of a SNP annotation tool BAD_Mutations that makes use of a likelihood ratio test based on alignment of all currently publicly available Angiosperm genomes.